Circuit QED lattices: towards quantum simulation with superconducting circuits 
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The Jaynes-Cummings model describes the couphng between photons and a single two-level atom 
in a simplified representation of light-matter interactions. In circuit QED, this model is implemented 

by combining microwave resonators and superconducting qubits on a microchip with unprecedented 
experimental control. Arranging qubits and resonators in the form of a lattice realizes a new kind 
of Hubbard model, the Jaynes-Cummings-Hubbard model, in which the elementary excitations are 
polariton quasi-particles. Due to the genuine openness of photonic systems, circuit QED lattices 
offer the possibility to study the intricate interplay of collective behavior, strong correlations and 
non-equilibrium physics. Thus, turning circuit QED into an architecture for quantum simulation, 
i.e., using a well-controlled system to mimic the intricate quantum behavior of another system too 
daunting for a theorist to tackle head-on, is an exciting idea which has served as theorists' playground 
for a while and is now also starting to catch on in experiments. This review gives a summary of the 
most recent theoretical proposals and experimental efforts in this context. 



I. INTRODUCTION 

Circuit QED describes the interaction between nonlin- 
ear superconducting circuits [1 3] and microwave pho- 
tons stored in resonators or propagating in transmission 
lines. To minimize Ohmic dissipation, the circuits used 
in this context arc made of superconducting materials. 
The currents and voltages in such a circuit are gov- 
erned by a quantum Hamiltonian and a corresponding 
Schrodingcr equation. Similar to atoms, superconduct- 
ing circuits have a ground state and low-lying excitations 
which form a discrete and anharmonic energy spectrum. 
In much of the research geared towards solid-state quan- 
tum information processing, the energy levels above the 
lowest two can be ignored and, in slight but common 
abuse of language, the circuit is referred to as a super- 
conducting qubit. Once coupled to a resonator, external 
microwave generators and detectors, the quantum state 
of the circuit can be manipulated and read out at will. 

Strikingly, such superconducting circuits exhibit quan- 
tum coherence on a macroscopic scale: the size of su- 
perconducting qubits and transmission line resonators 
ranges from hundreds of micrometers to millimeters. 
Nowadays, experiments with such circuits are routinely 
used in producing entangled states and in produc- 
ing effective photon-photon interactions, rendering them 
promising building blocks for larger quantum networks. 
This appeal of electric circuits, their modularity and rel- 
ative case of fabrication have greatly boosted the interest 
in the relatively young field of circuit QED [4-7]. 

The simplest circuit QED system consists of a sin- 
gle superconducting qubit coupled to the electromagnetic 
field of a single mode inside a microwave resonator (see 
Fig. 1, left panel). It represents an open-system real- 
ization of the well-known and exactly solvable Jaynes- 
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Cummings model [11] described by the Hamiltonian 

Hjc = oJrO^a + u)qa~^a~ + g{a^cr~ + aa~^). (1) 

Here, creates a photon with energy oJr in a particular 
resonator mode, are the raising and lowering opera- 
tors for the superconducting qubit with transition energy 
u)q, and g determines the strength of the qubit-photon 
coupling. Control and readout of the superconducting 
qubit can be modeled by incorporating driving and pho- 
ton leakage into an open-system version of the Jaynes- 
Cummings model (JCM). Such open-system formulations 
describe the reduced density matrix of the system as ob- 
tained from a Lindblad master equation or its unravelling 
in terms of quantum trajectories [12]. Natural extensions 
of the JCM including the Rabi model [13-15], the Dicke 
model [16, 17] and the Tavis-Cummings model [18-21] 
provide the framework for recent circuit QED experi- 
ments involving multiple qubits. At the time of writing, 
the state-of-the-art in circuit QED entails the implemen- 
tation of one, two and three-qubit gates on chips with 
up to four qubits [22, 23] and up to five resonators [23]. 
With coherence times of superconducting qubits now ap- 
proaching the 100 /iS mark [24, 25] , simplified versions of 
Shor's algorithm and quantum error correction have been 
demonstrated experimentally [22, 23]. 

The relatively slow rise of the number of qubits per chip 
-an increase from 1 to 4 in the last eight years (for com- 
parison: trapped ion experiments succeeded in demon- 
strating entanglement among 14 qubits in 2011 [26])- 
does not indicate limitations in fabrication. On the con- 
trary, fabrication of hundreds of qubits on a single chip 
is well within the realm of existing technology. The slow 
increase in qubit number rather points to the formidable 
challenges posed by quantum computation, in particular 
for long coherence times, very high gate fidelities, and 
individual addressability of qubits. From this viewpoint, 
building a circuit-based quantum simulator [27] rather 
than a quantum computer, might offer a less ambitious 
but, arguably, equally exciting goal. 
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FIG. 1. Circuit QED realization of the Jaynes-Cummings model and the Jaynes-Cummings lattice. Left: The device image 
(top) and circuit (bottom) show a transmission line resonator ("cavity") capacitively coupled to a nonlinear superconducting 
circuit ("qubit"). In the example shown here, the circuit is a transmon qubit [8, 9]. The transmon energy levels are tunable via 
magnetic flux provided by a flux-bias line ( "FBL" ) . "In" and "out" ports connect the system to microwave drive and detection 
circuitry. Device image and circuit from [10]. Right: Circuit QED realization of the Jaynes-Cummings lattice. Transmission 
line resonators, a, are coupled to each other via coupling capacitors, b, and thus form a regular lattice in which photons can 
hop from site to site. Photon-photon interactions are induced by the presence of superconducting circuits such as the transmon 
qubit, c. 



Indeed, well-controlled quantum systems with tunable 
Hamiltonian and many degrees of freedom might ul- 
timately provide invaluable clues about the nature of 
phases and dynamics of strongly correlated many-body 
systems for which the validity of various approxima- 
tion schemes not infrequently remains a matter of de- 
bate. Spectacular results have already been achieved 
with quantum simulators based on systems of ultracold 
atoms [28, 29] and trapped ions [30]. The idea that 
photon-based systems could likewise serve as quantum 
simulators goes back to the theoretical proposals by Hart- 
mann [31], Angelakis [32], and Greentree et al. [33]. 

Photons are an unconventional choice for a quantum 
simulator. As opposed to atoms, their number is not 
strictly conserved and their interaction is virtually neg- 
ligible for frequencies easily accessible in the lab. The 
Jaynes-Cummings model (1) indicates how to cure the 
latter ailment. In this model, the interaction between 
photons and matter leads to mediated photon-photon in- 
teraction. The circuit QED implementation of a photon- 
based quantum simulator thus consists of the step from 
a single Jaynes-Cummings system (see Fig. 1, left panel) 
to an array of coupled Jaynes-Cummings systems (see 
Fig. 1, right panel). The physics of the resulting Jaynes- 
Cummings Hubbard model (JCHM) 

HjcHM = J2 ^JC,j + J2iJalaj + h.c), (2) 

thus involves interaction on each site j described by the 
on-site Jaynes-Cummings term i?jc,j, and hopping of 
photons between nearest-neighbor sites with rate J. 

We will review both the evident similarities with the 
Bose-Hubbard model [34] as well as the important differ- 
ences, which could render circuit QED lattices a promis- 
ing quantum simulator for dissipative phase transitions 



[35]. Several excellent reviews of the field already exist 
[36-39] and are recommended to the interested reader. 
The present article aims to give a synthesis of the most 
important and most recent results of the field, with an 
emphasis on ideas advancing the circuit QED-based im- 
plementation of interacting photon lattices. The material 
we present is structured as follows. 

In section II, we introduce the generic model of a cir- 
cuit QED lattice and show how it reduces to the standard 
Jaynes-Cummings-Hubbard model. We discuss its valid- 
ity and the typical energy scales associated with its input 
parameters and show how the openness of the system is 
taken into account within a master equation formalism. 
In section III, we explain how effective photon-photon in- 
teractions arise in circuit QED due to light-matter cou- 
pling. We discuss the photon blockade effect, which is 
the basic underlying mechanism for the generation of 
strongly correlated photon states. In section IV, we sum- 
marize recent proposal and results for small, finite-size 
lattices. Finally, section V deals with large extended 
arrays and, in particular, focuses on the superfluid-to- 
Mott-insulator phase transition of polaritons. Our re- 
view closes in section VI with a summary and outlook on 
the future of circuit QED lattices. 



II. FROM CIRCUITS TO LATTICE 
HAMILTONIANS 

Living up to their name, circuit QED systems are con- 
veniently designed and described as electric circuits. The 
Hamiltonian of such circuits can be derived systemat- 
ically via circuit quantization [40-42]. We review the 
central steps and discuss their application to the primary 
building blocks of circuit QED lattices: superconducting 
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qubits and transmission line resonators. 



A. Circuit quantization 

Circuit quantization essentially follows the usual pro- 
cedure of canonical quantization, consisting of 

1. Formulation of the classical circuit Lagrangian, 

2. Legendre transform to the classical Hamiltonian, 

3. Quantization via p„ — >■ — «g§-, 

where the Xn and Pn denote generalized variables and cor- 
responding conjugate momenta. The majority of work is 
invested into the first step. To formulate a circuit La- 
grangian, we first select generalized coordinates {(j>n} for 
the circuit. The most convenient choice is to associate 
each node in the circuit (i.e., a conductor connecting cir- 
cuit elements to each other) with one generalized coordi- 
nate 0,1 (i) = /j*^ dt'Vn{t'). Here, Vn denotes the electric 
potential at node n. The resulting generalized coordi- 
nate is the time integral of the node voltage, starting at 
a fixed reference time tQ, and has dimensions of magnetic 
flux. For a circuit with N individual nodes, gauge free- 
dom allows us to eliminate one variable by letting the 
corresponding node act as ground, say = 0. 

Employing these coordinates and their associated ve- 
locities, we can specify all kinetic and potential energy 
contributions and express the Lagrangian in the usual 
form C{(f>m 4>n) = T — V . The common elementary build- 
ing blocks of superconducting circuits are capacitors, in- 
ductors, and Josephson junctions (JJ). For each such ele- 
ment, when connecting node n to the corresponding 
energy contributions are 



^capacitor — 2 ^'n) 
Vinductor — 2J^(^^+'^ ^ 



c 



L 



and 



-Ejcos \2'k( 



The special role of magnetic flux in superconductors pro- 
duces one additional rule: any superconducting loop (in- 
terrupted, at most, by Josephson junctions) is subject 
to fluxoid quantization and is influenced by any exter- 
nal magnetic flux <I>oxt threading the loop. The circuit 
responds to external flux with a loop current that ren- 
ders the total flux through the ring an integer multiple of 
the superconducting flux quantum $o = h/(2e). In the 
Lagrangian, this constraint is accounted for by selecting 
one element in each loop to act as the "closure branch," 
for which the variables enter as 







^)-*c 



(3) 



thus including the external flux. After Legendre transfor- 
mation and canonical quantization, one thus obtains the 



quantum Hamiltonian for the circuit. A word of caution 
is in order regarding "shortcuts" that claim to directly 
obtain the Hamiltonian: the Legendre transform, in gen- 
eral, requires the inversion of the capacitance matrix and 
is trivial only for the most simplistic circuits describable 
by a single generalized coordinate. 

Superconducting qubits. The circuits used as su- 
perconducting qubits [1-3] are prime examples for circuit 
quantization at work. The hamiltonians obtained using 
the above three steps allow for detailed predictions of 
energy spectra and other observables in good agreement 
with the host of experimental data on the various flavors 
of charge [9, 43-45], flux [46-48], and phase qubits [49]. 
In most cases, the "qubits" actually have excited states 
beyond the | ) and 1 1 ) states which are commonly ne- 
glected in the qubit hamiltonian 



Q 



E 



(4) 



Here, denote the Pauli raising and lowering operator 
for qubit n. Higher levels may be important to take into 
account for qubits with weak anharmonicity, such as the 
transmon qubit [9, 50]. In this case, black-box quantiza- 
tion introduced recently by Nigg et al. [42], provides an 
elegant alternative to conventional circuit quantization. 
In systems with closely spaced superconducting qubits, 
additional terms describing direct qubit-qubit interaction 
need to be included in Hq and could produce interesting 
"photon-less" lattices of superconducting qubits [1]. 

Transmission line resonators. Circuit quanti- 
zation as described above can also be applied to 
transmission line resonators (figure 2), and entail the 
usual lumped circuit discretization into LC elements 
and the subsequent continuum limit. Considering the 
associated equations of motion, one obtains the normal 
mode eigenvalue problem 
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with boundary conditions Tdx^tj.\x^ = ■^C't'^^'^a'Ut 
the left (— ) and right (+) endpoints of the resonator. 
The mode functions are normalized according to 



dx^lix) = l. (6) 



In these expressions, fj, enumerates normal modes, Cip 
are the capacitances at the left and right endpoints, and 
c and £ denote the capacitance and inductance per unit 
length. For each mode, ^fj.{x) gives the voltage at posi- 
tion X inside the resonator. The corresponding normal 
mode frequencies are obtained as the non-zero solu- 
tions to the transcendental equation 



tan a;,. 



1 



(7) 



where = Lx\/£cuj^ and Xt = C^/{cLx). = 
|a;+ — is the total length of the resonator. De- 
tailed derivations of these relations can be found in Refs. 
[51, 52]. 
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FIG. 2. Capacitively coupled transmission line resonator, a, Lumped element representation of resonator with inductance i and 
capacitance c per unit length. The resonator is connected to the remainder of the circuit on the left and right by capacitances 
C- and C+. b, Lowest three normal modes of a weakly coupled resonator and c, corresponding frequencies from graphical 
solution of equation (7) for x± ~ 0.07. d, Capactive three-way junctions with uniform coupling and experimental realization 
(lower panel, with permission, from [38]). 



Photon hopping. Arrays of resonators are formed 
by interconnecting resonators via coupling capacitors 
(figure 2d). As long as their capacitances are small com- 
pared to the total resonator capacitance, C± <^ cL, pho- 
tons can hop from one resonator to its nearest neighbors 
only. Focusing on a single mode in each resonator, the 
photon sublattice is thus described by 



Hn = 



J nn' ^n^n' 



(8) 



As usual, the operators aj^, a„ raise and lower the number 
of photons in resonator n with frequency ujr,n by 1 and 
obey the canonical commutation rules [a„,a3„] = 5nm for 
indistinguishable bosons. 

For full-wavelength mode and uniform hopping, Cc = 
C±, the hopping amplitude 



(9) 
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carries a positive sign [52]. Whenever half- wavelength 
modes are involved, the sign change of individual mode 
functions <i>'-"^ may induce negative hopping amplitudes 
and lead to modifications of the photonic spectrum if 
the lattice contains loops consisting of an odd number 
of resonators [51]. More advanced schemes with inter- 
mediate coupler circuits have been proposed for render- 
ing photon-hopping tunable in-situ and for inducing com- 
plex hopping phases associated with broken time-reversal 
symmetry [51, 53]. 

Qubit-photon interaction. The interaction be- 
tween qubits and resonator photons typically takes the 
form 



(10) 



assuming that each site of the circuit QED lattice con- 
sists of a simple Rabi or Jaynes-Cummings system, i.e., a 



single resonator coupled to a qubit located inside the res- 
onator. Generalizations of this, such as multiple qubits 
coupling to the same local resonator, are naturally con- 
ceivable and have been demonstrated experimentally. 
The specific origin and strength of this coupling differs 
according to qubit type and coupling scheme. For a 
simple superconducting qubit — namely, a single Joseph- 
son junction with capacitive shunt as used for transmon 
qubits — Devoret et al. [54] give a careful and instructive 
discussion of possible coupling types and their respective 
strengths. Capacitive coupling of the junction electric 
charge to the voltage field inside the resonator resembles 
the dipole coupling of an atom to the electromagnetic 
field. However, the small mode volume of chip-based 
transmission line resonators and the large effective dipole 
moments of superconducting qubits can render the effec- 
tive coupling strength g/tOr ~ 2l3{E,j /2EcY/'^y/a quite 
large. (Here, symbols denote a: fine structure constant, 
(3: capacitance ratio of order unity, Ej: Josephson en- 
ergy, Ec'- charging energy.) Values of g approaching the 
ultra-strong coupling regime (5 ~ cj^, Wg) are feasible by 
coupling the qubit directly to the current through the 
center-pin of a coplanar waveguide resonator [54-56] . 



B. Open-system description of circuit QED lattices 

Experimentally, circuit QED devices are operated at 
sub-Kelvin temperatures - not only to maintain com- 
monly used metals like aluminum and niobium deep in 
the superconducting phase but also to minimize spuri- 
ous population of the microwave resonators with ther- 
mal photons. In striking contrast to most other quantum 
simulators, circuit QED arrays realizing the JCHM have 
a rather trivial ground state, namely the vacuum state 
void of photons and with all qubits occupying the ground 
state. Here, a key difference from ultracold atoms arises 
from the fact that the grand-canonical ensemble does not 
exist for photons: coupling to a heat bath alone, i.e.. 
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Table 1 Qrcuit QE) parameters and typical values 



Model E^rameter m 


Typical values m/Zn 


ln-9tu tunability? 


Wr 


resonator frequency 


1... 15GHz 


Nfes- with SQUID terminated resonators [ 57] 




qubit frequency 


100MHz. ..15GHz[ 64] 


Nfes- e.g., via global or local magnetic flux [ 58] 


g 


photon-qubit coupling strength 


0... 400 MHz 


>fes[55,59] 


J 


photon hopping strength 


1...100MHz[ 60] 


■Vfes-with additional coupler circuits] 51 ,53] 


K 


photon escape rate at port 


10kHz... 80 MHz[60] 


>fes-with additional coupler circuits] 51 ,61] 


Tk 


intrinsic photon escape/dissipation rate 


> 5 kHz [62] 


No 


Ti 


qubit relaxation rate 


> 10 kHz [24] 


Nfes-e.g., via Purcell e ect [ 63] 


JM 


qubit dephasing rate 


> 50 kHz [24] 


>fes- by inducing extra noise ^^^^^ 



working in the canonical ensemble, already introduces ad- 
ditional photons into the system. (An alternative word- 
ing describing the same situation is that the chemical 
potential for photons always vanishes, see [65].) 

A significant host of theoretical work has been based 
on the assumption that a chemical potential might still 
be engineered if the time scale for reaching a quasi- 
equilibrium with given photon number is much shorter 
than the time scale at which photons escape from the 
system. The resulting simplicity of equilibrium physics 
at T = is highly convenient for the theorist. Whether 
this regime is accessible in the experiment, however, re- 
mains an open question. One way or another, the full 
theoretical description of circuit QED arrays will need 
to account for the essential openness of this quantum 
system: the array is driven out of equilibrium with mi- 
crowave generators to introduce photons, and photons 
leaking out of the array reveal crucial information about 
the quantum state inside. 

Provided coupling to the environment is sufficiently 
weak and power spectra of associated noises show little 
variations at the frequencies of interest, Lindblad mas- 
ter equations and their unraveling in terms of quantum 
trajectories constitute a robust framework for studying 
dynamics as well as steady-state physics of driven arrays 
[66, 67]. A fairly generic master equation for the reduced 
density matrix p describing a circuit QED array then 
takes the form 

dtp = -z[i?JCHM + -ffdrivo, P] + 71 X! '^^'^n]P 

n 

+ 7^ ^ I?[<]p + l^Y. ^t""]'^ + '^11 'DMp. 

n n n'Gports 

(11) 

Here, TJdrive = ^(ame'"''* -I- aj^e"*""**) denotes classi- 
cal drives with frequency Ud and drive strength ^ which 
act on all or a subset of all resonator sites. As usual, 
the damping operators are defined as 'D\L\p — LpL^ — 
^{L^Lp + pL^L). The terms on the right-hand side of 
equation (11) describe unitary evolution of the system 
under its Hamiltonian i/jcHM + ^^drive, spontaneous re- 



laxation and pure dephasing of qubits with rates 71 and 
7ip as well as photon loss from the array due to unwanted 
dissipative channels (rate 7^) as well as intended loss of 
photons at ports of the array (rate k) which are then 
directed to detectors. For ports realized by capacitive 
coupling of an array resonator to an external transmis- 
sion line, the intended loss rate is given by the expression 
K = AZqC^uj]^ [60], where Z^^flfc is the characteristic 
impedance and Co the output capacitance at the port. 

For reference, we summarize the most important pa- 
rameters for circuit QED arrays along with typical ranges 
for their values in table II A. Entries do not claim funda- 
mental limitations but rather reflect theoretical proposals 
and current experimental capabilities. 

III. PHOTON BLOCKADE 

Effective photon-photon interaction can be realized 
with photons in a cavity, which strongly couple to a com- 
ponent of matter, e.g., via dipole coupling to a two-level 
system (qubit) as described by the JCM in Eq. (1). The 
eigenstates of the JCM are superpositions of a Fock state 
with n photons and the qubit occupying the ground state 
(g), or (n — 1) photons with the qubit occupying the ex- 
cited state (e), yielding 

= sin(?„|n ,5) + cos6'„|(n — 1) , e) , 
\n~) = cos.en\n,g) - fiin6n\{n - 1) ,e) . (12) 

Here, we have introduced tan(0„) = 2g^/n/ {5-\-2xn) with 
Xn = \/ g^n + <5^/4 and the detuning parameter 5 — ujr~ 
ujq. The corresponding eigenvalues are 

= - 5/2 ± Xn ■ (13) 

The ground-state \0,g) with energy eq = is the only 
state that does not hybridize. In Eq. (13) upper and 
lower polariton energies are separated by the Rabi split- 
ting 2xn and, on resonance, have the linewidth Se = 
(71 + 27^ + 7,)/2. 

Polaritons can thus be interpreted as photons dressed 
by the qubit through coherent emission and absorption 
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6.9 

tJd/27r [GHz] 

FIG. 3. Supersplitting of the vacuum Rabi resonance beyond 
linear response as presented in [68]. The bottom panel shows 
the vacuum Rabi splitting as measured with low drive power. 
For sufficiently high drive power, the single Lorentzian splits 
into a doublet of peaks as shown in the top panel, (with 
permission, from [68]). 



processes, i.e., Rabi oscillations. The dressed photon 
then inherits the anharmonicity of the matter component 
leading to a nonlinearity in the spectrum of the combined 
light-matter system, which is given by the energy differ- 
ence 



u = i4 



l) - (4 



(14) 



The parameter U quantifies a spectral shift which de- 
scribes the energy cost of adding a second photon to the 
cavity versus adding the first. It is sometimes also called 
particle hole gap or effective Hubbard-U, resembling the 
solid-state analogue (we emphasize that this analogy does 
not hold for higher excitations n > 2, where the nonlin- 
earities of the Hubbard model and the JCM significantly 
differ). The spectral shift in Eq. (14) is largest for zero 
qubit-cavity detuning (i5 = 0), where U = g{2 — \/2). 
It becomes vanishingly small for large detunings, i.e., 
U = 0{g^/6^). In this dispersive regime photons and 
qubits barely interact with each other. The nonlinearity 
in Eq. (14) has been observed using direct pump-probe 
spectroscopy [68, 69] as well as collapse and revival ex- 
periments [70]. 

One of the most spectacular consequences of the non- 
linearity in Eq. (14) is the photon blockade effect [71, 72]. 
When a weak coherent drive is resonant with the lowest 
transition in the spectrum [tud ~ — eo) the cavity will 
be populated with a single photon, but the entrance of a 
second photon is effectively blocked. This photon block- 
ade is the fundamental mechanism behind the realization 
of strongly correlated photon states. Its observation in 
various cavity QED systems [73, 74] has triggered an im- 
mense interest in using coupled light-matter systems for 
quantum simulation. We next summarize the signatures 
of photon blockade in single-cavity circuit QED systems. 



Using a coherent microwave drive, signatures of pho- 
ton blockade in circuit QED systems have been observed 
in measurements of heterodyne transmission, resonance 
fiuorescence spectra, as well as the photon statistics in 
linear response. A necessary requirement for the clear 
manifestation of photon blockade is that the spectral 
shift in Eq. (14) is larger than the drive strength and 
the linewidth of the polaritons {U ^ £,,6e) so that multi- 
photon processes are suppressed. 

For a weakly driven Jaynes-Cummings system, the 
heterodyne transmission T ^ ] (a) ] exhibits Lorentzian 
peaks when the drive frequency matches with lower or 
upper polariton frequency. Beyond linear response, the 
Lorentzian splits and develops a central dip, i.e, an an- 
tiresonant lineshape termed supersplitting in Ref. [68] 
(see Fig. 3). This phenomenon can be simply under- 
stood in the truncated Hilbert space of only two states in 
the Jaynes-Cummings ladder, namely the vacuum state 
and the upper or lower single-polariton state. Within 
this two-level system, the heterodyne transmission signal 
is equivalent to the polarization of the dressed two-level 
system and can be calculated exactly [68]. For strong 
driving, the transmission vanishes exactly at the center 
of the linear-response transmission peak. This somewhat 
surprising result originates from the saturation of the 
two-level system in an equal mixture of its two states, 
rendering the average polarization and thus the hetero- 
dyne transmission signal effectively zero. 

Inside a nonlinear cavity, the qubit first gets dressed 
by the quantum cavity field (leading to the effective two- 
level system described above) . Adding an additional clas- 
sical drive thus corresponds to a dressing of dressed states 
[71]. Another manifestation of this dressing is the ap- 
pearance of a triplet peak structure in resonance fluores- 
cence spectra, known as the MoUow triplet [75-78]. 

One of the most striking consequences of photon block- 
ade can be seen in the statistics of the emitted light. In 
the blockade regime, no more than one photon is present 
inside the cavity at a time. Only if this photon is emitted, 
e.g., leaving through the leaky mirrors, another photon is 
allowed to enter. When the time the photon spends inside 
the cavity is much larger than the decay time, i.e., un- 
der strong coupling conditions, the photon decay leads to 
an anti-bunched train of photons leaking out of the cav- 
ity. A nonlinear cavity operating in the photon blockade 
regime can thus be regarded as a single photon turnstile 
device. The strong suppression of the zero-time delay 
second order correlation function [781, 



g(2)(0)<l with g('^>{t) 



(at a) 2 



(15) 



is a tell-tale sign of the photon antibunching. 

An alternative scheme for observing photon blockade 
in circuit QED consists of exciting the system with a 
broadband microwave signal [59] instead of the conven- 
tionally used coherent tone. By varying the spectral 
bandwidth of such an incoherent drive, the transmitted 
power displays several steps similar to what is known 
as the current staircase for electron transport through 
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FIG. 4. Upper panel: (a) Coulomb blockade in a quan- 
tum dot with discrete A^-electron resonances separated by the 
charging energy. Electron transport through the dot is con- 
trolled through the gate voltage, Vg, and source-drain voltage, 
Vad- (b) Photon blockade in a nonlinear cavity with discrete 
A'^-polariton resonances separated by the Jaynes-Cummings 
nonlinearity. Photon transmission through the cavity is con- 
trolled by the center frequency, lOd and the bandwidth of the 
incident photons, a [59]. Lower panel: Calculated number 
of photons inside the cavity n for typical experimental pa- 
rameters [59]. The inset shows a measured photon staircase, 
i.e, the normalized transmission T, for fixed ujd as a function 
of bandwidth a analogous to what is known as the Coulomb 
staircase in mesoscopic electron transport, (with permission, 
from [59]). 



Coulomb blockaded quantum dots (see Fig. 4). Here, the 
spectral bandwidth of the impingent signal plays a role 
analogous to bias voltage in mesoscopic electron trans- 
port [79, 80]. 

Strong nonlinearities at the single photon level can also 
be realized in setups with a continuum of electromag- 
netic modes - for example, by embedding quantum two 
or few level systems directly into waveguides or trans- 
mission lines [81-83]. In this case, embedded system can 
act as a scattering center for an incoming train of pho- 
ton wave packets. Due to the interaction, photons are 
either transmitted or reflected. Here, the coupling of a 
few discrete states to a continuum of propagating photon 
modes (rather than coupling to a discrete mode localized 
in a cavity) constitutes a major difference to conven- 
tional cavity QED setups. Typical nonlinear quantum 
optical phenomena such as electromagnetically induced 
transparency (EIT) and photon blockade may be real- 



ized in such systems without the need of high-Q cavities 
[84-86]. Experimental progress realizing this new strong 
light-matter coupling regime has also been made with 
devices coupling metallic or semiconductor nanowires to 
quantum dots [87, 88] , diamond nanowires to NV centers 
[89] , and cold trapped atoms to the continuum of modes 
inside a hollow fibre [90]. 

Within circuit QED, experiments with superconduct- 
ing qubits coupled to open one-dimensional transmission 
lines have led to observations of stronly suppressed trans- 
mission and the MoUow triplet in the RF spectrum of 
scattered photons [77, 91]. Strong quantum correlations 
resulting in the anti-bunching of reflected photons have 
also been observed using transmon qubits [92]. These 
first experiments might set the stage for an interesting 
new research area, in which quantum impurity phenom- 
ena such as many-body bound states [93, 94], or the 
Kondo effect [95] can be studied by optical means. 



IV. SMALL LATTICE SYSTEMS 

Small lattices of Jaynes-Cummings systems can be con- 
sidered as fundamental building blocks in quantum net- 
works with applications in quantum computation and 
communication. Circuit QED technology is a promising 
approach in this context due to the high level of coher- 
ent control, low-noise properties and potential scalability 
on a microchip. While single-qubit manipulations are 
routinely realized using microwave drives, a major chal- 
lenge is the controlled generation of entanglement among 
many qubits. So far this has been demonstrated in circuit 
QED with several qubits inside a single cavity using var- 
ious qubit-qubit coupling mechanisms [58, 96-99]. Two 
and three-qubit gates have been built [22, 23, 58, 98, 100] 
and entanglement between three qubits and the genera- 
tion of GHZ states with high fidelity has been achieved 
[101, 102]. A remaining important task is the demonstra- 
tion of non-local entangled states between two or more 
qubits in distant cavities using circuit QED arrays with 
a few resonators [32, 103-105]. 

Progress in this direction has been achieved by en- 
gineering a controlled coupling between two and three 
empty resonators using qubits as inter-cavity couplers 
[106-108]. As a key step towards the realization of 
circuit QED based quantum simulators experimental- 
ists have recently reported the first successful fabrica- 
tion of a twelve-resonator array arranged in a Kagome 
geometry [60]. By careful design of coplanar waveguide 
and characterization of multiple samples via transmis- 
sion measurents, disorder in resonator frequencies was 
successfully pushed down to the level of one part in 
10^. From a theoretical point of view, small lattice sys- 
tems are numerically tractable and exact solutions for 
dynamics and steady states have been obtained in vari- 
ous contexts[31, 32, 50, 52, 109-116]. These theoretical 
studies may serve as a testbed for future experiments. 

What is the fate of the photon blockade in small cou- 
pled systems, where the nonlinearity of the spectrum 
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may strongly depend on the size of the array A'"s and 
the strength of the hopping parameter J ? The level 
scheme for the lowest excitations of a Jaynes-Cummings 
dimer (consisting of two coupled Jaynes-Cummings sites) 
as well as the band structure of large lattices are shown 
in Fig. 5. Choosing all qubits to be resonant with the 
symmetric photon state of the dimer or the bottom of 
the photon band in large arrays ensures that the qubit 
subspace remains strongly coupled to the cavity array as 
the hopping strength J is increased. 

For weak hopping (J ^ g) the Jaynes-Cummings array 
exhibits a photon blockade similar to the one expected for 
a single Jaynes-Cummings site. In this regime, Jaynes- 
Cummings arrays may be utilized as quantum simulators 
of various spin lattice systems, where virtual photon hop- 
ping mediates XY or ZZ couplings between the qubits 
[32, 109]. An interesting application is the simulation of 
the phase diagram of high-spin Heisenberg models [117] 
which are very difhcult to calculate numerically. Recent 
proposals suggest that JC lattices could also be used 
to simulate effective spin models with topological non- 
trivial states, e.g., Majorana-like bound states [118-120]. 
Even though these states would not be topologically pro- 
tected, circuit QED may offer a unique detection scheme 
for a proof of principle of their existence. 

For strong hopping (J ^ y), photon blockade sensi- 
tively depends on the number of cavities in the array. 
To leading order in 5/ J, the nonlinearity of a ID chain 
with periodic boundary conditions and qubits in reso- 
nance with the bottom of the photon band, is given by 
[121] 

U[N,) = 2g(l - Vl - l/(2iV,)) + Oigyj) . (16) 

Consequently, the nonlinearity decreases with the num- 
ber of sites in the array. For instance, a dimer produces 
C/(2) « g{2 -V5)+ Oig^/J) while U{1) = g{2 - \/2) for 
a Jaynes-Cummings site. 

The condition g ^ Se, J for strong coupling, however, 
is not always necessary for the manifestation of photon 
blockade in coupled Jaynes-Cummings arrays [122, 123]. 
If a Jaynes-Cummings dimer is driven asymmetrically 
such that only one of the two cavities is pumped, anti- 
bunching can occur even under weak-coupling conditions 
( g Se,J). The reason for this unconventional ap- 
pearance of blockade is a subtle quantum interference 
of photons hopping between the cavities. Observing 
anti-bunching in this regime is more difficult because 
the second-order correlation function defined in Eq. (15) 
oscillates rapidly on a time scale set by the hopping 
strength J. Detection would thus require fast photon 
detectors [123]. The first single photon detectors in the 
microwave regime were tested only recently [124-126]. In 
the opposite regime of ultra-strong coupling where the 
qubit-cavity coupling becomes comparable to the cavity 
and qubit frequencies {g 0;^,^^), photon blockade is 
substantially modified due to the presence of strong two- 
photon decay processes [127]. 

The competition between photon blockade and pho- 
ton hopping in small lattices has also been studied in the 
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FIG. 5. Upper panel: Energy level diagram for the 1- 
excitation manifold of the JCHM for the dimer {Na = 2) and a 
large array (A^s = 00) [121]. The qubit (Q) energy is kept res- 
onant with the symmetric cavity state (S) of the JC dimer and 
the bottom of the photon band of the array {cuq = oJr — J)- The 
drive frequency ujd matches the bottom of the lower polariton 
band iPk {oJd ~ uJr — J—g)- Lower panel: Second-order coher- 
ence (?'^'(0) for zero time delay of the dimer (exact numerics) 
and an infinite array (mean-field approximation) as a function 
of the hopping strength J/g [121]. The drive strength is fixed 
to //g = 0.001 (black), 0.005 (red), 0.01 (green), 0.02 (blue), 
(with permission, from [121].) 



context of photonic Josephson-like effects [112, 114]. Mo- 
tivated by analogous work on weakly interacting atomic 
[128-130] and polaritonic systems [131, 132], photons in a 
Jaynes-Cummings dimer have been proposed to undergo 
a sharp self-trapping transition from a regime where an 
initial photon population imbalance between the two res- 
onators undergoes coherent inter-cavity oscillations (de- 
localized) to a regime where it becomes self-trapped (lo- 
calized) as the photon-qubit interaction is increased [114]. 
A coherently driven Jaynes-Cummings array with three 
sites has also been proposed to act as a Josephson in- 
terferometer [112]: when the outer sites are driven pho- 
ton blockade manifests in a suppression of Josephson os- 
cillations in the center cavity, leading to characteristic 
changes in the statistics of emitted photons. 

The impenetrable photon regime in a ID chain repre- 
sents an extreme case of a strongly-correlated many-body 
state of photons. Absorption spectra are predicted to ex- 
hibit unique signatures of photon fermionization into a 
Tonks- Girardeau gas under non-equilibrium conditions, 
for small arrays with as few as only five sites [113]. In 
related work, the Tonks-Girardeau regime has been sug- 
gested to be accessible by engineering dissipative pro- 
cesses [115]. Utilizing such reservoir engineering for driv- 
ing the Jaynes-Cummings array into a specific steady 
state may also enable the formation of nonequilibrium 
photon condensates inside arrays with only weak inter- 
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FIG. 6. Mean-field phase diagram for tlie JCHM display- 
ing the lowest three Mott lobes with polariton numbers 
= 1,2,3 [154]. Dotted lines represent the critical hopping 
strengths Jc/ff for which chemical potential fi and detuning 5 
are chosen such thatl particle-hole symmetry is obeyed. Finite 
detuning > decreases the critical hopping strength Jc/g 
for iV > 1. For the lowest Mott lobe (iV = 1), Jc/g steadily 
increases when increasing the detuning S. (with permission, 
from [154].) 

actions [133]. 

V. LARGE LATTICE SYSTEMS 

Large lattice systems allow to study the interplay of 
collective behaviour in extended systems and strong cor- 
relations due to local interactions. For example, inter- 
acting bosons on an infinite lattice as described by the 
Bose- Hubbard model (BHM), undergo a superfluid-to 
Mott insulator (SF-Ml) phase transition at zero tem- 
perature [34]. Motivated by the experimental realiza- 
tion of this transition with ultra-cold atoms in an op- 
tical lattice [135], the study of strongly correlated and 
condensed phases of ultra-cold atoms has become an ex- 
tremely active and versatile research field over the last 
decade [28, 134]. 

The recent realization of strong light-matter interac- 
tions in cavity and circuit QED systems has likewise 
triggered an immense interest in realizing such strongly 
correlated phases and condensates in photonic systems. 
Bose-Einstein condensates (BEC) of weakly interacting 
polaritons have already been demonstrated experimen- 
tally with exciton-polaritons inside a quantum well when 
coupled to a photonic crystal cavity [136, 137]. Spec- 
tacular experimental advances have shown the existence 
of superfluidity [138, 139], quantized vortices [140], and 
quantum solitons [141] in these systems. While the equi- 
librium versus non-equilibrium nature of these new quan- 
tum condensates remains a matter of debate [132, 142- 
144], a key experimental challenge now is to reach the 
limit of strong correlations, e.g., a possible SF-MI tran- 
sition of polaritons [31-33]. 

Circuit QED may offer a well-suited approach to reach 



this regime since nonlinearities on the few photon level 
are routinely realized in single cavity systems. The de- 
sign of large lattice systems as described by the JCHM 
in Eq. (2) seems also feasible with current state of the 
art technology [38]. Below we summarize the extensive 
theoretical work in this context. 



A. Quasi-equilibrium 

For an infinite array of coupled Jaynes-Cummings 
sites, each comprised of a single photonic mode coupled 
to a two-level system, the Jaynes-Cummings-Hubbard 
Model (JCHM) predicts a SF-MI phase transition of 
polaritons [32, 33]. The competition between strong 
photon-qubit coupling, giving rise to an effective repul- 
sive nonlinearity (localization), and the photon hopping 
between cavities (delocalization) leads to an equilibrium 
quantum phase diagram featuring Mott lobes reminis- 
cent of those for ultracold atoms in optical lattices as 
described by the BHM [34]. 

Similar to the BHM, the JCHM obeys a global U{1) 
symmetry associated with the conservation of the total 
number of polaritons N. In order to calculate the JCHM 
phase diagram it is thus convenient to work in the grand- 
canonical formalism and to introduce a chemical poten- 
tial /i for polaritons according to Hjchm ^^jchm— M-^- 
Results for the quantum phase diagram can then be de- 
rived from a simple mean-field theory by decoupling the 
hopping term via [33] 

J2al<^j ~ z{a^ij + a^p* -\^p\^) . (17) 

Here we assume an infinite system with lattice transla- 
tional invariance and denote the resulting order param- 
eter by V' = (a) ■ The coordination number z gives the 
number of nearest neighbor cavities per site. This decou- 
pling scheme maps the original lattice model to an effec- 
tive local on-site Hamiltonian. The phase diagram shown 
in Fig. 6 is then obtained by minimizing the ground-state 
energy of the local mean-field Hamiltonian with respect 
to ip. 

Inside the Mott lobes, on-site repulsion dominates over 
photon hopping and local fluctuations of the photon num- 
ber are strongly suppressed corresponding to a localized 
Mott state for photons and polaritons (since the qubits 
are trivially localized). Outside the Mott lobes, photons 
are delocalized over the lattice and form a weakly inter- 
acting polariton superfluid. The quantum phase transi- 
tion between these two phases takes place at some critical 
value of the hoping strength Jc/g (at fixed detuning 5), 
where the Mott lobes terminate. 

Formally, this mean-field approximation only becomes 
exact in the limit of large coordination numbers, z oo 
[145]. In practice, the mean-field diagram shown in Fig. 6 
is in qualitative agreement with exact numerics down to 
z — A, e.g., for a two-dimensional square lattice. For 
one-dimensional lattices the transition is of BKT type 
and cannot be described by mean- field theory [146] . The 
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precise location and shape of the phase boundary can 
be determined numerically by using various methods in- 
cluding DMRG calculations in ID [146, 147], variational 
cluster approximations [148] and Quantum Monte Carlo 
simulations [149-152]. Modified models capturing next- 
nearest-neighbor, diagonal and long-range photon hop- 
ping in ID and 2D also lead to substantial changes in 
the Mott lobe sizes [152]. 

Important differences between the BHM and the 
JCHM arise from the composite nature of polariton 
quasi-particles. Depending on the value of the detuning 
parameter S, the nature of the polaritons can be photon- 
like {6 < 0) or qubit-like (5 > 0) . At the same time, finite 
detuning S ^0 weakens the nonlinearity of the spectrum 
and thus reduces the size of the Mott lobes - with one 
exception: the size of the lowest Mott lobe (iV = 1) in- 
creases for positive detuning. Only in this case the nature 
of all polaritons in the lattice become qubit-like and are 
thus trivially localized. 

An analytic strong-coupling theory for the phase dia- 
gram and the elementary excitations in the Mott phase 
of the JCHM has been derived based on a linked-cluster 
expansion (LCA) [145]. The random-phase approxima- 
tion (RPA) of this LCA reproduces the mean-field phase 
diagram and predicts the doubling of gapped excitation 
modes in the Mott phase as compared to the BHM [153] . 
The effect of quantum fluctuations, i.e., corrections to 
the RPA, have also been calculated analytically using the 
same method [153] and are in rough agreement with exact 
results [148, 149, 152]. In the superfluid phase, a slave- 
boson theory [154] and an effective action approach [155] 
as well as numerical methods [149] predict the existence 
of a gapless linear Goldstone mode and gapped amplitude 
or Higgs modes in close analogy to the BHM [156]. While 
the Goldstone mode was observed in trapped weakly in- 
teracting atom [157, 158] and polariton superfluids [139], 
a measurement of the Higgs mode has only recently been 
reported with strongly correlated atoms trapped in an 
optical lattice [159]. 

The RPA [153] and the slave-boson approach [154] as 
well as simple field-theoretic considerations [160] predict 
that the JCHM changes the universality class at the tip 
of the Mott lobes, where it belongs to the 3D XY uni- 
versality class. This has recently been confirmed using 
large-scale Monte-Carlo simulations [151]. 

The SF-MI transition of polaritons described above 
gets substantially modified in the ultra-strong coupling 
regime. In this case, counter-rotating terms like aV+ 
and aa^ have to be included in Eq. (1). These terms 
break the U{1) symmetry of the JCHM. They can be ne- 
glected as long as g <^ oJr,ojq, but need to be included 
for ultra-strong couplings with g ^ u)^.,ujq. In the pres- 
ence of these counter-rotating terms, the JCHM maps 
to the transverse field Ising model, which features a dis- 
crete Z2 parity symmetry-breaking transition [161]. The 
same phase transition could also be realized in a chain of 
Josephson junctions with direct electric dipole coupling. 
Here, circuit QED can be used for high-resolution spec- 
troscopy by attaching transmission line resonators to the 



end of the chain [162] . 



B. Non-equilibrium 

The assumption of an effective chemical potential in 
the previous section reqiiires that the system equilibrates 
on a time-scale much shorter than the lifetime of po- 
lariton quasi-particles. While this appears applicable 
to some experimental setups with large coherence times 
such as phonon- polaritons in trapped ion systems [163], 
it is yet unclear whether the same can be achieved within 
the circuit QED architecture. Taking into account the 
basic non-equilibrium nature of most quantum optics ex- 
periments, in particular the presence of drives and dissi- 
pation, is one of the biggest challenges in the field. New 
tools and methods are needed to investigate the inter- 
play between strong correlations and collective behavior 
under non-equilibrium conditions. 

For the Bose- Hubbard model, quench dynamics was 
studied to assess the effects of dissipation on the SF-MI 
phase transition [164] . Initializing the system with a 
suitable excitation pulse in a Mott-like state, the sys- 
tem response turns out to be highly sensitive to the 
properties of the underlying equilibrium ground-state. 
While the absolute value of the order parameter decays as 
]V'(i)l ~ e~'*'~* and thus contains no valuable information, 
the rescaled coherence parameter V'(t) — \ip{t)\/ \Jn{t) re- 
veals the equilibrium SF-MI phase transition unambigu- 
ously. In particular, in the long-time limit one obtains 
i){t — >■ 00) ~ for J < Jc, i.e., in the Mott phase, and 
^j{t 00) > for J > Jc, i.e., in the superfluid phase 
( Jc denotes the critical hopping strength of the quantum 
phase transition in equilibrium at the tip of the lobe). 
The distinction between Mott and superfluid phases be- 
comes sharp only when 7^ — >• 0, but is smoothened and 
washed out when 7^ is non-zero. In fact, dissipation 
tends to favor the localized regime [164, 165]. Conse- 
quently, quench dynamics allow for the observation of 
clear and sharp signatures of an equilibrium phase tran- 
sition even in open dissipative systems. 

Qualitatively new effects can be expected when dissi- 
pation is compensated by an external pump, which drives 
the system into a non-equilibrium steady state, rather 
than the ground state. Large arrays with up to 16 cavi- 
ties, each of them driven by a coherent source, have been 
treated in [166] using a time evolving block decimation 
algorithm for the BHM. A similar algorithm has recently 
also been applied to the JCHM [167]. This numerical 
technique allows to treat relatively large system sizes but 
is limited to a regime with small inter-site entanglement, 
i.e., for weak photon hopping. In Ref. [166] it was shown 
that when a weak drive is applied to each cavity with 
opposite phases on neighboring sites, photons crystallize 
into a dimerized steady state with positive correlations 
in neighbouring cavities, but anti-correlations for larger 
separations. 

Steady states in a driven, infinite lattice have been 
considered in Ref. [121] using a generalization of the de- 
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coupling mean-field approximation in Eq. (17) to den- 
sity matrices. This yields a nonlinear, but closed sys- 
tem of equations for the time evolution of local operator 
expectation values. For weak hopping, one finds that 
an infinite JC lattice remains photon blockaded. For 
strong hopping, photon blockade breaks down and key 
signatures like anti-resonant lineshapes in homodyne de- 
tection, Mollow triplets in RF spectra, and antibunch- 
ing of photons are destroyed [121]. Note that this is in 
agreement with the perturbative result in Eq. (16), which 
predicts that the nonlinearity vanishes in the thermody- 
namic limit {Ns — > oo). 

The crossover from weak to strong hopping is smooth 
as long as the array is pumped at the bottom of the 
lower polariton band [121]. Thus, no traces of the equi- 
librium quantum phase transition remain, which is to be 
expected since the external drive explicitly breaks the 
U{1) symmetry of the Hamiltonian. However, higher 
pump frequencies may cause tunneling induced bistabil- 
ities [121]. Contrary to the usual bistable behavior in 
quantum optics known for a single cavity [168-171] and 
for the Dicke model [16], here, the bistability develops 
at fixed pump strength when increasing the hopping rate 
J between cavities on a large lattice and might thus be 
interpreted as a first-order non-equilibrium phase transi- 
tion from a localized to a delocalized state. 

So far, the JCHM was discussed in the context of a pos- 
sible superfluid-Mott insulator transition of polaritons. 
However, the JCHM also represents the multi-mode ver- 
sion of the Dicke model, which is commonly used to de- 
scribe weakly interacting polariton EEC's in solid-state 
micro-cavities [142, 144]. In circuit QED, the single-mode 
Dicke model [16] can be realized in a chain of qubits cou- 
pled near resonantly to a single mode of a transmission 
line resonator. It is a current matter of debate whether 
one can realize a superradiance transition with such a 
device, if operated in the ultra-strong coupling regime 
[14, 17]. 

Jaynes-Cummings lattices also show rich and versatile 

physics in the context of artificial gauge fields for pho- 
tons [51, 52, 172-174]. By now, multiple mechanisms 
have been proposed to break time-reversal symmetry in 
circuit QED lattices [51, 52, 174]. For instance, time- 
reversal symmetry may be broken by using additional 



passive superconducting circuits as inter-cavity couplers 
instead of employing direct capacitive coupling. This 
design gives rise to complex phases in the photon hop- 
ping terms, which can be used to study further inter- 
esting strongly correlated phenomena including topolog- 
ical phases and chiral edge modes [175] or the fractional 
quantum Hall effect [176, 177]. First experiments in this 
direction are currently underway [174]. 

VI. CONCLUSIONS 

With the remarkable level of control over small circuit 
QED systems and the first experimental realizations of 
small circuit QED lattices at hand, the field is rapidly de- 
veloping and attracting physicists with a variety of back- 
grounds ranging from quantum optics and mesoscopic 
physics to condensed matter physics and strongly corre- 
lated systems. Combining circuit QED with other archi- 
tectures in so-called hybrid systems has already begun 
and will further continue to drive the field [178]. 

Prom the theoretical perspective it is important to 
identify model systems simple enough for obtaining re- 
liable, quantitative predictions. Such models will serve 
as a testbed for the characterization and optimization of 
the first experimental implementations of circuit QED 
lattices. At the same time, new mathematical tools and 
approximation schemes need to bo developed in order to 
treat non-equilibrium physics, collective behavior of large 
ensembles and strong correlations on an equal footing. 
Of course, the ultimate goal of a quantum simulator is to 
invert the roles of theory and experiment: experimental 
results for large systems will provide crucial information 
and guidance for the development of robust theoretical 
models of interacting quantum systems with many de- 
grees of freedom in and out of equilibrium. 
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